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Abstract 

Accurate analytical expressions for the state densities of liquid 4 He droplets are 
derived, incorporating the ripplon and phonon degrees of freedom. The microcanon- 
ical temperature and the ripplon angular momentum level density are also evalu- 
ated. The approach is based on inversions and systematic expansions of canonical 
thermodynamic properties. 



1 Introduction 

Important dynamical processes in finite systems such as nuclei, polyatomic molecules, 
nanoclusters, atomic clouds, droplets frequently turn out to be statistical in nature: 
evaporation/fragmentation, radiation, emission of electrons, equilibration between inter- 
nal degrees of freedom or between host and solvent molecules. When such a system 
is thermally isolated, e.g. when flying in a beam or suspended in a trap, the proper 
statistical-mechanics treatment is that of the microcanonical ensemble where the energy 
E is fixed and not the temperature of an external heat bath. The density of states, or 
level density, p(E)dE represents the number of quantum states between energy E and 
E + dE. For separable degrees of freedom this is number of normal mode combination 
such that their energies add up to a total internal energy lying in this interval. The func- 
tion plays a crucial role in the thermal description of microcanonical systems. For low 
excitation energies p(E) can be represented by a sum of delta functions, corresponding 
to excitations of a only a few of the individual modes, but for even moderate excitation 
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energies the density of these delta functions becomes so large that is is well described as 
a continuous function of energy. In this situation it is most convenient to use a density 
smoothed over the discreteness of the energy levels. In addition to energy systems de- 
scribed with the microcanonical ensemble have a conserved total angular momentum, so 
the correspondingly resolved density of states, p(E, J), is often of relevance. 

Free liquid helium nanodroplets [H [2] represent an interesting system for a statistical 
treatment. One reason is that helium is the only element which cannot be described in 
terms of classical dynamics for any internal degrees of freedom under the experimental 
conditions used to study the droplets. This makes the system interesting in its own right. 
For example, "magic number" maxima in the size distributions of small 4 He clusters have 
been shown to correlate with the ability of the cluster to accommodate elementary ex- 
citation modes |3j. A second reason is the use of the droplets as micro-cryostats used 
to investigate other clusters and molecules. Evaporative cooling generates internal ener- 
gies corresponding to temperatures of ~ 0.4 K and is used to thermalize impurities to 
this otherwise unreachable temperature for gas phase molecule and cluster beams. Under- 
standing these processes requires accurate density-of-states expressions for the elementary 
excitation spectrum. 

The calculation of level densities requires that the excitation spectrum is known. At 
low temperatures, the relevant normal modes of 4 HeAr clusters within the liquid drop model 
are ripplons which are quantized capillary surface waves, and phonons which are quantized 
bulk compression waves. For large droplets these modes are separable to a good approxi- 
mation [1] , a fact that greatly facilitates a statistical analysis of the excitation spectrum. 
For a spherical droplet, both ripplon and phonon modes possess well-defined eigenvalue 
spectra characterized by angular momentum for ripplons, and angular momentum and 
mode index for phonons. A calculation of the total density of states requires enumeration 
of all possible normal mode combinations, with individual energies and angular momenta 
adding up to a given total E and J. 

A leading-order evaluation of p{E) for ripplons was carried out by Brink and Stringari 
[5]. Subsequently, Lehmann [6] presented a comprehensive discussion of the densities of 
states for ripplons and phonons computed by direct numerical counting, and showed that 
the resultant plots of the logarithm of the level densities could be well parameterized by 
polynomial fits. These fits were then used to calculate other thermodynamic functions 
and to analyze droplet cooling with angular momentum conservation constraints [S] • 

In this paper, we show that accurate density of states functions can be obtained by 
analytic evaluation. This is appealing in its own right, as the calculations take advantage 
of several elegant and generally useful tools from the literature. In addition, having ana- 
lytic expressions for various types of elementary excitations provides a systematic method 
for treating situations where several types of normal modes are excited simultaneously, or 
when the spectrum of elementary excitations is modified. 

The plan of the paper is as follows. In Section [2] we calculate the ripplon density of 
states as a function of energy, Section [3] considers its angular momentum dependence, Sec- 
tion H] is devoted to phonon excitations, Section [5] to the angular momentum of phonons, 
and Section O to the total p(E) function. Section [7] comments on the similarity between 
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the spectra considered here and those of multielectron bubbles in bulk liquid helium, and 
presents a summary. 



2 Ripplon density of states 

As mentioned above, ripplons are quantized waves on the droplet surface. For a spherical 
liquid drop, the elementary excitation spectrum is given by [9] 

e t = tw oy /e(e-i)(e + 2). (1) 

Here i > 2 is the angular momentum quantum number of the wave and 

where a t is the coefficient of surface tension, D the mass density, R the droplet radius, m a 
the atomic mass, and TV the number of atoms in the droplet. If the parameters of bulk 
liquid helium are used, we have Huq ~ ^3.8/\ZiV^K in temperature units [6]. Below, the 
ripplon energy will be expressed in dimensionless units scaled to the quantity Tiloq. Each 
mode has a degeneracy of {21 + 1). 



Canonical approximation 

A first approximation to the level density can be derived in the canonical ensemble picture, 
where it is assumed that the system possesses a definite temperature T, and the system's 
internal energy is associated with its most probable value. The energy density of states 
of a finite system is then given by [TTJJ, [TTJ [T2] 

p(E) = &S (3) 

'-2ir{dE/d0) 



where (3 = (k B T) 1 and 

S = pE + lnZ (4) 

is the entropy; Z is the canonical partition function. In the following we will use units 
where ks — 1- The square root appearing in the equation involves the heat capacity 
and appears because the canonical entropy includes an approximately gaussian integral 
over the thermally populated states with a width given by the heat capacity and the 
temperature, see, e.g., Ref. [T2] . 

Since the ripplon elementary excitations are bosons we have 

lnZ = - ]T(2^+l)hi(l-e- fe ). (5) 

t=2 
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The canonical thermal energy of the ripplon ensemble is 



E = —d(lnZ)/d/3. (6) 

To leading order, we can replace the sum in Eq. (jSj) by an integral from zero to 
infinity, and approximate the energy eigenvalues (fTT) by Be & £ 3 / 2 . The integral then 
straightforwardly evaluates to 

\nZ = r (I) c (1) P~ 4/3 = 1.685/?- 4 / 3 , (7) 



,3/ ' \3, 

and from Eq. ([6]) the (dimensionless) energy is 

E = 2.247/T 7/3 . (8) 

Assembling everything into Eq.([3]) and expressing the answer in terms of the energy, 
we find 

p r i P (E) « 0.31LET 5/7 exp(2.476£ 4 / 7 ), (9) 
which is the same answer as in Ref. [5]. 

Micro canonical ensemble 

The above calculation can be improved in two places. One obvious refinement is to 
evaluate the sum (jSJ) with greater care and to use more precise eigenvalues. A deeper 
conceptual question is how to compute thermodynamic quantities for a finite isolated 
system for which the total internal energy is a conserved quantity and not an expectation 
value and the use of a "temperature" must be carefully defined. A thorough discussion was 
given by Andersen et al. in Ref. [12] with the conclusion that the convenient canonical 
formalism may be retained, but the canonical expression for the energy must be 
corrected as follows: 

E = -d(\nZ)/dp-p-\ (10) 

Here E is the fixed excitation energy of the system and j3 is understood as the "micro- 
canonical temperature" defined as 

f3 = d[lnp{E)}/dE. (11) 

The procedure taken is as follows. First, the sum in Eq. <$±j§ is calculated using the 
first three terms of the Euler-Maclaurin summation formula [13J. With the form of the 
spectrum given in Eq.([T]), this formula becomes 

-\nZ = y"(2£+l)ln(l-e- fe )^ + -ln(l-e- /3e2 ) (12) 
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The upper limit, £ max , has been set to infinity as before. The actual value is on the order of 
P-max ~ 2ttR/ X m in ~ 2nR/ (2d), where A is the wavelength and d is the interatomic distance 
[I]. In the liquid drop approximation (R = N 1 ^ 3 d/2) one then has £ max ~ nN 1 ^ 3 /2. In 
view of Eqs. (00E]) this yields a size-independent ripplon Debye temperature of e max ~ 7.5 
K. Using this value to estimate the error in In Z, the leftover terms are found to be on the 
order of (Pe max /4: — 7£ max /6) exp(—/3e max ). For T = IK this is a relative contribution to 
ln(Z) of less than 10~ 2 /iV 1//3 which will be ignored. 

A tedious calculation of Eq. (fl2l) involving expansions of exponentials in powers of (3 
results in [TJ] 

349 7 

In Z = 1.685/T 4/3 + 0.639/T 2/3 - — + - \n(2V2(3) + ... (13) 

96 3 

The first term coincides with Eq. ([7]), and the rest are finite-size and spectral correc- 
tions. Note that all the numerical coefficients derive from explicit expressions involving 
special functions. The expansion Eq. lfTB"]) has been checked against a numerical sum. The 
comparison is shown in Fig. HJfor Helmholtz' free energy, F = —T\n(Z). Already at tem- 
peratures where T is equal to the lowest excitation energy £ 2 = hcu ^/6, the free energy 
is well represented by the above expression. At higher energies the agreement improves 
monotonically. 

Knowing the partition function, we can now use Eq. fflUj) to determine the relation 
between the microcanonical energy and temperature: 

E = 2.247/T 7/3 + 0.426/T 5/3 - — /T 1 . (14) 

3 

Again, the first term reproduces Eq. (jEJ). 

In order to proceed with the calculation of the entropy, the heat capacity and the 
level density in Eq.flBJ), we need to invert the relation (TI4~I) which expresses E(f3) to get 
(3(E). This is done by the iterative method of successive approximations. The result is 
an expansion for in powers of E~ 2 /"' ' , 

ft- 1 = 0.7069£ 3/7 - 0.07239£ 1/7 + 0.7212£~ 1/7 + (15) 

where the coefficients are calculated from those entering Eq.(ll4p. Now the prefactor and 
the exponent in Eq.(j3]) can be evaluated, using Eqs. fj4^. fTl3|) .fri4l). and [3(E), finally 
yielding 

pn p (E) = 0.205E- 12/r exp (2.476£ 4/7 + 0.507£ 2/7 ) . (16) 

Let us emphasize again that all the numerical coefficients encode analytical expressions. 

Eq. (fl6l) . which is the main result of this section, may be compared with the canonical 
approximation (Q, an exact numerical count carried out with the help of the Beyer- 
Swinehart algorithm [T5], and the form written down in Ref. [B] as an empirical fit to 
the numerical count in the interval E=50-2500. Fig. [2] shows such a comparison, and 
demonstrates that the analytical expression gives an excellent representation of the exact 
result [16]. 
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3 Ripplon angular momentum density 



The next step is to generalize the ripplon state density to a function which is not only 
energy- but also angular momentum-resolved. This problem has been comprehensively 
studied in nuclear physics [10, E] ■ One way of visualizing the net angular momentum of 
a large distribution of excitations with varying (£, £ z ) is as the result of random angular 
momentum coupling, in which case the central limit theorem applies and one expects to 
find a normal distribution. Indeed, the above references show that p(E, J) is essentially 
a product of p(E) and a Gaussian factor: 

^■ J ^^2^-^- (17) 

It is permissible here to replace J(J + 1) by (J + \) 2 . 

It is still necessary to establish the variance a 2 . An elegant way to do this to leading 
order by means of an extended grand canonical distribution is described in Bethe's review 
[TU] . where the method is applied to a system of non- interacting fermions in a spheri- 
cal potential box. Here we follow the same procedure for a system of bosonic ripplon 
excitations. 

The idea is first to evaluate the projection M of the net angular momentum J of the 
droplet in terms of the contributions of individual normal modes at temperature T. The 
fact that M is a conserved quantity is accounted for by a separate Lagrange multiplier, 
or "chemical potential" 7, such that 

7 = -OS/dM (18) 
(S is the entropy). We can calculate M directly by summing over all modes: 

^ * m r , „ r m • dm , , 

1=2 m=—£ o — £ 

(in reduced energy units). Expanding the integrand to first order in 7 [TO] , we find 

M - l r (I) < (!) ^ 1/3 < < 20 > 

and Eq. (jSJ) allows us to express the result in terms of the droplet energy. To leading 
order we have: 

7 = l.77QME~ s/7 . (21) 
Now we can use Eq. (|T8|) with Eq. ()2ip to obtain the entropy variation: 

S(E, M) = S{E, 0) - M 2 /(2a 2 ) (22) 

with 

(2(T 2 ) _1 = 0.888E" 8/7 . (23) 
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The second term in Eq. (j22p leads to a normal distribution in M. The distribution in 
J can be shown to have the same variance [TOl [TT] . Therefore Eqs. (1171) and (]2"3"]) define 

Prip(E, J). 

The numerical evaluation of the rotational density of states in Ref. [6] led to essentially 
the same form of the state density function, with the factor corresponding to (2<r 2 ) -1 fitted 
as 0.868E~ 8 / 7 + 0.964.E -13 / 7 , which affirms the analytical result (1231) : the factors deviate 
by less than 2% for E =100-2500. 

A shorter estimate of the variance is illustrative. The number of quanta in one mode 
(£, m) is on the order of T /ei for levels up to et ~ T and zero for higher quantum energies. 
The total number of excited quanta is then 

n « ( 2i + ^)T/£ 3/2 w 4T 4/3 , (24) 

where the sum was approximated by an integral and T is written in terms of the Uq 
unit. With the energy-temperature relation (fT4"]) we get the leading order value for en- 
ergy per quantum (e) = E/n = 2.247T/4, and from this an average of {£) = (e) 2//3 = 
T 2 / 3 (2.247/4) 2 / 3 . The standard deviation a of £ is then, according to the 'random walk' 
argument used above, a = \/n(£). Inserting the calculated (£) and expressing the result 
in terms of the total energy, one has 

(2a 2 ) _1 = 2- 1 / 3 (2.247)~ 4/21 E^ 7 = 0.68£- 8 / 7 , (25) 

in surprisingly sensible agreement with the above result. 

One may seek to describe the angular momentum distribution in the language of 
a rotational energy and a moment of inertia /, associating [TT] the exponential in Eq. 
(I17p with a Boltzmann factor involving /3h 2 J(J + 1)/ (2/), i.e., I = h 2 [3a 2 . Using the 
canonical-ensemble results, Eqs. OH]) and (123]) . we can express the "ripplon moment of 
inertia" in terms of the ripplon excitation energy (in reduced units) : 

/ = 0.797£ 5/7 . (26) 

4 Phonon density of states 

Surface ripplons are the lowest-temperature droplet excitations; bulk phonons appear 
next. These are compression sound waves which arise as solutions of the wave equation 
within the volume of the spherical drop. As such, their energies are given by 

e n/ = huk nt £ (27) 

where u is the speed of sound and the wave number k n ^ is determined by the boundary 
condition at the surface. If the Dirichlet boundary condition is adopted [HE], then /c n ^ = 
a n> e/R, where a n ^ is the nth root of the je spherical Bessel function. For a free surface, a 
more appropriate boundary condition is the Neumann one, in which case k n ^ = a' n£ /R, 
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with a' ni the root of the Bessel function derivative, j[. The energy scale is set by the 
longest wave length, i.e., k ~ tt/R, so we can express phonon energies in units of 

i = hwr/R, (28) 

which works out to e = (25 

.5JV 1 / 3 ^K in temperature units if the speed of sound in bulk 
4 He is used [B] . The leading-order behavior of the phonon state density can be determined 
in a straightforward way by invoking the standard expression for the Debye heat capacity 
(per unit volume) of bulk phonons: 

27r 2 k 4 

Cbulk = IbW^^' (29) 

Multiplying this by 47ri? 3 /3 and using the fact that (in the canonical framework) 
C = dE/dT and {k B T)~" = dS/dE, we can use integrations to express S in terms of E. 
Then, from Eq.([BD we find that to first order, In p P h(E) ~ S(E) = 3A1E 3 ^. This matches 
the leading term of the fit to a direct numerical count in Ref. [6] which is 3.331E 13 / 4 . 

The Debye temperature for phonons in liquid 4 He is ~ 25 K [21 J, corresponding to a 
total phonon thermal energy (from Eq l3Tl) of ~ lOOOiV K. We can therefore use the low 
temperature approximation throughout. 

The prospect of refining the calculation by analytically evaluating a statistical sum over 
the precise spectrum (1271) may seem bleak, as the Bessel function roots which "contribute 
in an essential manner are just the zeros for which the usual formulae (like McMahon's 
expansion) are bad approximations" jTTj. However, rescue comes from an elegant math- 
ematical framework known as the Weyl expansion [T71 [HI IS] ■ It provides a systematic 
expression for the smoothed density of eigenmodes in a finite cavity in terms of volume, 
surface, and curvature terms. As described in the above references, this is a very general 
theory, valid for both scalar and vector wave equations, and applicable to a wide variety 
of physical phenomena. 

Ref. [20J applied this formalism to the specific heat of metal nanoparticles. The finite- 
size correction to the specific heat (|29|) derived there is immediately usable for our droplet 
problem: 

r-r <\ 9C(3) k BT\iklT 

c ~ Cbulk H —J^li + gmW (30) 

The + sign applies to the Neumann and the - sign to the Dirichlet boundary conditions. 
Although we focus on the Neumann condition, the Dirichlet case will be included for 
completeness. 

We now follow almost the same sequence as in the bulk limit described above: Eq. 
( 1301) is multiplied by the droplet volume and integrated once to obtain (with E and T in 
units of e) 

E(T) = ^T 4 (i) C(3)tt 2 T 3 + ^T 2 , (31) 
45 9 
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and a second time to obtain S(T) as S = / 7 (C ' /T')dT' '. The first function is inverted by 
iteration to yield 

T(E) = 0.391£ 1/4 (+) 0.069 - O.OOG.Er 174 . (32) 

[Calculating T(E) instead of [3(E) is more convenient in this case.] Eq. (j3J) is then 
used to obtain the density of states. The calculation is done to the first three orders 
in E, in correspondence to the three terms in the heat capacity expansion (|30|) . The 
microcanonical correction (fTUj) in the present case turns out to contribute only in the 
next order of smallness. The result of the calculation is as follows: 

p ph (E) = AE~ 5/S exp ^3.409£ 3/4 (-) 0.908£ 1/2 + 0.482£ 1/4 ) (33) 

Once again, the + sign is for the Neumann boundary condition on the phonon wave at the 
droplet surface and (-) for the Dirichlet condition. Using the bulk canonical heat capacity 
in Eq. fl30l) gives a pre-exponential factor of A = 0.26. 

FigE] compares the exact Beyer-Swinehart count for the phonon spectrum with the 
full Eq.fl33l and with the level density based on the bulk Debye heat capacity, Eq. fl29l) . 
i.e., where only the first term in the exponent is present. Fig.fll]) shows a more detailed 
comparison of Eq. fl33|) and the exact-count phonon level density. We find good agreement 
between analytical expression and the exact computation, although not as good as for the 
ripplon case. 

The estimate of the prefactor A in Eq. (13"3"|) cannot be expected to be correct because it 
does not include higher-order expansion terms in the exponent that would yield corrections 
of the same order. A comparison with the numerical count suggests a correction in the 
form of a factor exp(— 0.62E 02 ). Although this correction is larger than the error found 
for the ripplon level density, it is nevertheless still relatively small. An effective value of 
A pa 0.05 can be used for energies below 400. 

5 Phonon angular momentum density 

A computation of the angular momentum resolved phonon level density suffers from the 
difficulties with expressing the lowest Bessel function eigenvalues with a simple functional 
form. In contrast to p p h(E) there is, to our knowledge, no solution in the literature 
for this problem. As will be clear from the results presented in section [6] below, the 
contribution to the level density from the phonons is minor compared to that of the 
ripplons, and the required precision in the calculation of the angular specified phonon 
contribution is therefore correspondingly smaller. In this section we will make an order 
of magnitude estimate, based on the leading order term of McMahon's expansion of the 
roots of the Bessel functions |13j . For the Neumann boundary condition the roots are 
(n + £/2 — 3/4)7r pa (n + £/2)ir. With the phonon energy scale used, Eq. ([28l . the quantum 
energies are thus n + £/2. When states with energies up to T are averaged, the linear 
dependence of the quantum energy on £ gives an average value of (£) ~ T. Since also the 
n-dependence is linear the constant of proportionality is on the order of unity. The total 
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number of states below energy T is on the order of T 3 . Combining these estimates give, 
using the same type of 'random walk' estimate as Eq. 0251) for the ripplons, that 

~ ^ = (!0 5/4 ^ 5/4 - 10 °£- 5/4 - ( 34 ) 

The ratio of the cr's for the phonons and ripplons (here denoted a r ip) with the leading 
order terms in the caloric curves Eqs. fll4|3ip and the proper energy scaling is 

^ ~ 0.002 (T[K]) 7/e iV 1/6 . (35) 
(jj.jp 

This is small compared to unity up to extremely large droplet sizes. The conclusion 
that the width of the phonon angular momentum distribution can be ignored holds very 
well, even if the estimate of the width should be incorrect by as much as an order of 
magnitude. 



6 Combined level density 

A helium droplet may have both ripplon and phonon oscillations excited at the same time 
(and, at higher temperatures, rotons as well pQ). The coupling between these normal 
modes is weak at bulk liquid surfaces [22] and in large droplets [I], thus their energy 
contents may remain independently defined for some length of time, and the individual 
state densities will then come from Eqs. ([TBI) and ([3"31 . The question of equilibration 
dynamics of excitations in superfluid droplets and the relevant time scales is a very in- 
teresting one, and has not yet been addressed in detail. Below, we discuss an estimate 
of state densities in circumstances when the ripplon and phonon excitations do achieve 
statistical equilibrium. 

In principle, the level density of combined excitations can be calculated by direct 
summation, as described in Section [2J This would be a very involved procedure, because 
the ripplon and phonon quantum energies have different dispersion relations and scale 
differently with size. Alternatively, one can calculate the level density as a convolution. 
Also in this task does one benefit from formulating the general problem in terms of the 
microcanonical temperature. The convolution to be performed is 

p(E) = / p rvp (E - e)p ph (e)de. (36) 

For not extremely large droplets the largest part of the excitation energy resides in the 
ripplons. Indeed, the ratio between the energies of the ripplon and phonon subsystems is, 
canonically: 

|^ « 6.8 x 10~ 3 iV 1/3 (T[K}f 3 . (37) 

(Temperature expressed in Kelvins.) It is clear that for temperatures under 1 K (i.e., 
those which lie safely below the Debye cut-off values specified above and below the onset 
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of roton modes) and droplets of up to several tens of thousands of atoms in size, the phonon 
energy contents is a fraction of the ripplon energy. Under these conditions one can treat 
the ripplon degrees of freedom as a heat bath and calculate the phonon contribution with 
an expansion of the integrand of Eq. fl36|) around some small phonon energy. We will use 
the simplest choice of zero phonon energy, and to increase the precision we expand the 
logarithm of the level density. Thus 



P{E) 



p r ip(E) exp 



d\n[p rip (E)] 1 2 d 2 \n[p ri p(E)] 
£ dE + 2 £ dE 2 



p ph {e)de. (38) 



The upper limit of the integral can be replaced by infinity without serious loss of precision 
because the integrand peaks well below this value. We recognize the first derivative in 
the exponential as the microcanonical temperature 1/T of the ripplon system at energy 
E, see Eq. (fTTj) . and therefore have 



p(E) = p r i P (E) e~ £/T pp h (£)exp 
Jo 



1 2 d 2 \n[p rip (E)} 
2 £ dE 2 



de. 



(39) 



The second exponential in the integrand can be expanded, with the integral of the first 
term yielding the phonon canonical partition function at T, Zp^T): 



p(E) = p np (E) lz ph (T) + £ 



e- £/T Pph (e) 



1 2 d 2 \n[p ri p(E)} 
2 £ dE 2 



de 



(40) 



To leading order and ignoring the difference between the canonical and microcanonical 
temperatures, this simplifies to 



p(E) = p ri p{E)Z ph {T) 1 - 



a 



ph 



/? 2 

^ph 



2 Crip 



2C r ipT 2 



+ .. 



(41) 



Hence the ratio of the term which is second order in e to the zero order term in Eq. ((40 
is approximately 



a 



ph 



2C„ 



+ 



E ph 



rip 



2C T ; L pT 



6 x 10- 3 (T[K]) 5/3 N 1/3 + 4 x 1(T 6 (T[K]) 1A/3 N 



4/3 



(42) 



For not excessively large or warm droplets we can leave out the correction terms and thus 
have 

p(E) = Prip (E)Z ph (T), (43) 

where Z p h(T) as stated above is the phonon canonical partition function at the micro- 
canonical ripplon temperature corresponding to the ripplon energy E. 

The exponential part of the phonon canonical partition function can be calculated, 
e.g., by integration of the standard relation Eq.([6]) with the caloric curve in Eq. (1311) . This 
procedure does not determine the integration constant which translates into a multiplica- 
tive constant on the total level density, Eq. fj43|) . This constant, c, is approximately the 
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product of the pre-exponential from Eg. (1551) and the prefactor that appears in Eq.flSJ (i.e., 
the value given by a saddlepoint expansion of the phonon level density in the calculation 
of the canonical particion function). The result is 



where C p h is again the phonon heat capacity. The leading order expressions for the phonon 
parameters C p h(T), E p h(T) give, taking into account the different scaling of energies for 
phonons and ripplons, the total level density 



p(E) = pHp(E) ■ 0.526iV 1/6 exp (o.04713A^- 1/2 T 3 + 0.01317iV 1/3 T 2 + 0.1634iV 1/6 T) (45) 



with the equation given in ripplon energy units and T = T(E) given by Eq. (fl5|) . The 
constant of 0.05 for the phonon level density pre-exponential, mentioned at the end of 
Sec. HI was used here also. 

This result is compared with numerical convolutions for N = 10 3 in Figj5j The 
numerical convolution was also calculated for iV = 10 4 with a similar result. 

One remark about Eqs. (14311451) is in place: These equations should only be used for 
calculations of microcanonical properties. For the calculation of canonical properties one 
should use the product partition function, Z riWiV h = Z rip Z p h. A naive application of 
Eq. ( 1451) in a calculation of the partition function of the combined ripplon-phonon system 
will give a divergent result at all temperatures. The origin of this divergence is the 
breakdown at high excitation energies of the approximations leading to the equation. 

7 Conclusions 

We have presented an analytical evaluation of the statistical density of states functions 
of the elementary excitations (surface ripplons and volume phonons) of isolated liquid- 
drop helium nanoclusters. These functions are expressed in terms of microcanonically 
conserved quantities: energy and angular momentum. The obtained formulas accurately 
match numerically computed curves as the energy level densities vary over ~ 150 — 300 
orders of magnitude. 

Other interesting helium systems to which the results may be applicable include 
micron-sized superffuid fog [23] and multielectron bubbles in liquid helium. The lat- 
ter are spherical voids inside bulk He, with a thin shell of electrons lining the inner wall 
(see, e.g., Refs. [2HE2] an d references therein). They can undergo small-amplitude shape 
oscillations, i.e., ripplons, whose frequency under zero external applied pressure has the 
form ujf oc (£ 2 — 1)(£ — 2) , which for large I approaches the same form as the droplet 
ripplon dispersion, Eq. ([I]). This implies that the statistical mechanics of these bub- 
bles should be similar to that of nanodroplets. One distinction is that the bubble are 
submerged into a bulk helium thermal bath, therefore for them the canonical ensemble 
treatment is rigorously correct and not just a convenient approximation. 

Finally, it should be pointed out that the results obtained in the present paper have a 
universal form and are expressed in terms of dimensionless scaled energies, therefore they 
are generally applicable to the statistics of droplets of various substances besides helium. 
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Figure 1: The (negative of) the ripplon free energies, calculated with Eq. (fl3l) (dotted line) 
and the summation in Eq. ([5]) (full line) which is exact apart from setting the upper limit 
to infinity. The temperature corresponding to the energy of the lowest excitation, £ 2 , is 
indicated. 
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Figure 2: Comparison of the ripplon level densities calculated according to Eq. (fT6!) (full 
line) and Ref.glj (dashed for E < 2500, dotted line for E > 2500). The fit in Ref . [6J was 
limited to energies between 50 and 2500, in the reduced units used here and is calculated 
as the derivative of the numerical fit to the integrated level density. The expressions have 
been divided by the exact Beyer- Swinehart result, causing the oscillatory behavior at low 
energy, and the curves plotted are the logarithms of these ratios. The curve of Ref. [5 J 
(not shown) is around 3. 
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Figure 3: Phonon level densities calculated according to the exact Beyer-Swinehart count 
(open circles), Eq. (l3"3l) (full line), and the level density derived from the bulk Debye heat 
capacity, i.e. corresponding to Eq.( l33l) but including only the first term in the exponential 
(dashed line). 
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Figure 4: A comparison of Eq. (l33l) and the exact-count phonon level density, showing 
essentially the relative difference in the entropy of the phonon system in the two calcula- 
tions. 
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Figure 5: The convoluted level densities for phonons and ripplons for droplet size 10 3 . 
The numerical convolution is the full line, and the approximate result in Eq.( H5|) the, 
hardly discernible, dashed line. The level densities for ripplons alone (dotted line) is 
given for reference. The arrow indicates the energy content of the ripplon excitations at 
a temperature of 1 K. 
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